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Single crystals of (Ndi_xCex) 2 Fei 4 B are grown out of Fe-(Nd,Ce) flux. Chemical and structural 
analysis of the crystals indicates that (Ndi_xCex) 2 Fei 4 B forms a solid solution until at least x = 0.38 
with a Vegard-like variation of the lattice constants with x. Refinements of single crystal neutron 
diffraction data indicate that Ce has a slight site preference (7:3) for the Ag rare earth site over the 
4/ site. Magnetization measurements show that for x = 0.38 the saturation magnetization at 400 
K, a temperature important to applications, falls from 29.8 for the parent Nd 2 Fei 4 B to 27.6 /is/f.u., 
the anisotropy field decreases from 5.5 T to 4.7 T, and the Curie temperature decreases from 586 to 
543 K. First principles calculations carried out within density functional theory are used to explain 
the decrease in magnetic properties due to Ce substitution. Though the presence of the lower-cost 
and more abundant Ce slightly affects these important magnetic characteristics, this decrease is not 
large enough to affect a multitude of applications. Ce-substituted Nd 2 Fei 4 B is therefore a potential 
high-performance permanent magnet material with substantially reduced Nd content. 

PACS numbers: 75.50.Cc, 75.50.Vv, 75.50.Ww, 71.20.Eh 


I. INTRODUCTION 

The reduction of critical materials present in consumer 
and industrial products in general, and in Rare Earth 
(RE) -based permanent magnets in particular, is a press¬ 
ing concern due to increasing demand and a decline 
and/or uncertainty in the supply chain of these materials. 
Dysprosium, judged to be the most critical of these ele¬ 
ments [1], is substituted onto the Nd site in Nd 2 Fei 4 B 
(2-14-1) permanent magnets to increase the value of 
their coercive field and maximum operating temperatures 
[2, 3]. Depending on the application, the amount of Dy 
added can be large, in the range of 1.4-8.7 wt% [1]. The 
highest temperature applications, those that rely on the 
greatest quantities of dysprosium, are required by the 
rapidly emerging technologies of wind turbines, magnet¬ 
ically levitated transport, and traction motors for hybrid 
and electric cars [1]. 

Though Dy is indeed considered the most critical el¬ 
ement, by many criteria Nd is not that far behind on 
the list [1]. It is therefore imperative that new and in¬ 
novative solutions are applied to the synthesis of perma¬ 
nent magnets that contain more abundant elements. To 
this end, recent work by Pathak et al. [4] has shown 
that partial substitution of Ce for Nd and Co for Fe 
{i.e. Ndi, 6 Ceo. 4 Fei 2 Co 2 B) results in a permanent mag¬ 
net with properties superior to those of Dy-substituted 
2-14-1 magnets for temperatures above 450 K. Since 1) 
Ce is by far the most abundant of the rare earth elements 
and 2) these new magnets simultaneously eliminate the 
need for Dy and reduce the amount of Nd used, they 
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represent an exciting opportunity for the development of 
lower-cost high performance permanent magnets. 

The crystal structure of Nd 2 Fei 4 B was established as 
tetragonal (space group P42/mnm, No. 136) as early 
as 1984 by the authors of Refs. [5] and [6]. There are 
two separate and inequivalent RE sites; one is a Ag site 
and the other a 4/ site. Also present are six inequiv¬ 
alent Fe sites and one unique B site. Past investiga¬ 
tions as to the effects of chemical substitution of dif¬ 
ferent RE elements onto the Nd site are numerous, see 
for example Ref. [7] and the references therein. How¬ 
ever, detailed, quantitative experimental studies on the 
effects of Ce substitution of varying amounts in pure sin¬ 
gle crystal samples have not previously been performed. 
Abache and Osterreicher [8] studied aligned powders to 
determine the crystal structure, magnetic anisotropy, and 
spin reorientation temperature in (Ndi_xCex) 2 Fei 4 B for 
X = 0.25 and x = 1 along with a multitude of other 
compositional variations of Nd 2 Fei 4 B achieved by vary¬ 
ing the degree and type of substitutions on the RE and 
transition metal (TM) sites. However, no detailed study 
was made concerning the magnetic properties as a func¬ 
tion of Ce concentration in this work. They also do not 
make any claims whatsoever as to the stability of the 
(Ndi_xCex) 2 Fei 4 B compound for 0.25 < a; < 1. More 
recently, density functional theory calculations of Alam 
et al. [9] actually predict a phase segregation for a; > 0.3 
in (Ndi_xCex) 2 Fei 4 B. Pathak et al. [4] have experimen¬ 
tally observed phase segregation under the conditions re¬ 
quired for the fabrication of melt-spun ribbons. 

Information also of interest in the (Ndi_xCex) 2 Fei 4 B 
system is the site preference of the Ce as there are two 
structurally distinct RE sites. Previous theoretical pre¬ 
dictions about Ce site preference exist [8, 9], but di¬ 
rect experimental evidence to determine whether Ce in 
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TABLE I: Molar Ratios for Sample Fluxes and Compositions 
of Single Crystals as determined from EDS. 

Flux Composition 

(molar percent) Composition of Composition of 

Nd Ce Fe 2-14-1 phase 1-4-4 phase 

51 0 43.3 5.8 Nd2Fei4“B NdFe4^^B4 

43.3 7.7 43.3 5.8 (Ndo.9i(8)Ceo.o89(7))2^®i4^^^ Ndo.904(9)Ceo.o96(i)Fe4^^B4 
32.9 18.1 43.3 5.8 (Ndo.78(4)Ceo.22(i))2l^®i4^^B Ndo.75(9)Ceo.25(3)Fe4^^B4 

25.5 25.5 43.3 5.8 (Ndo.62(6)Ceo.38(4))2Fei4^^B Ndo.6i(3)Ceo.39(2)Fe4^^B4 

18.1 32.9 43.3 5.8 “ Ndo.49(2)Ceo.5i(2)Fe4“B4 

“2-14-1 phase did not form large crystals 


(Ndi_xCex) 2 Fei 4 B is preferentially associated with the 
larger Ag site or the smaller 4/ site has previously been 
lacking. 

The (Ndi_xCex) 2 Fei 4 B system has recently come un¬ 
der much attention. Experimental investigations by 
Pathak et al. discovered an anomaly evident in the plot 
of the c lattice parameter vs. Ce concentration in melt- 
spun (Ndi_xCex) 2 Fei 4 B ribbons for x « 0.2 [4]. Sim¬ 
ilar anomalies were evident in plots of coercivity {He), 
maximum energy product {BH-eaax), and remnant mag¬ 
netization {Br) as functions of Ce concentration, along 
with the aforementioned phase segregation. The same 
investigation yielded an extraordinary result: when Co 
is substituted for Fe by two atoms per formula unit 
in (Ndi_xCex) 2 Fei 4 _yCoyB where x = 0.20, the Curie 
temperature (Tc) increased by 150 K while coercivity 
dropped only 22% as compared to an equivalent sample 
with no Co doping, i.e. x = 0.2, j/ = 0. The anisotropy 
field {Ha), saturation moment {Mg), and Br remained 
more or less unchanged compared to the x = 0.2, y = 2 
case. The possible applications of such excellent Dy-free 
magnets spawns an exciting line of research focused on 
linking the microscopic magnetic properties in pure crys¬ 
talline samples with the macroscopic magnetic properties 
of bulk permanent magnets. 

In the current work, we have employed single crys¬ 
tal growth techniques to synthesize single crystals of 
Ce-doped 2-14-1 compounds. We show here observa¬ 
tions contrary to those of the melt-spun ribbons of Ref. 
[4] for our samples grown under the slow crystalliz¬ 
ing conditions of flux growth; we are able to stabilize 
uniform single crystals of (Ndi_xCex) 2 Fei 4 B for values 
0 < X < 0.38. Furthermore, we have used these single 
crystal specimens to experimentally determine that Ce 
has a slight site preference in the 2-14-1 structure for 
the larger RE Ag site [8, 9]. Next, we have also mea¬ 
sured important magnetic properties such as Mg, Tq, 
Ha etc. for (Ndi_xCex) 2 Fei 4 B for the composition range 
0 < X < 0.38. Finally, we conclude this report with first 
principles density functional theory calculations which 
help explain the causes of the changes in Mg, Tc, and 
Ha with Ce substitution. 


II. METHODS 

Single crystals of both Nd 2 Fei 4 B and NdFe 4 B 4 type 
materials were grown from Nd-Fe flux following the tech¬ 
niques originally reported by Canfield et al. [10-12] and 
further refined by Saparov [13]. The starting materials 
were cuttings of high purity metals: Ce (Ames Labora¬ 
tory, 99.99%), Nd (Ames Laboratory, 99.99%), Fe (Alfa 
Aesar, 99.98%), Co (Alfa Aesar, 99.95%) and isotopically 
pure ^^B (^^B, ORNL, 99.99-1-%). The use of ^^B was 
necessary for neutron diffraction experiments as ^*^B is 
exceptionally good at neutron capture. Appropriate sto¬ 
ichiometries of these elements (see Table I) were loaded 
into Ta crucibles (1.25 cm diameter, 7 cm length) and 
sealed under 0.5 atm Ar using an arc-melter. A Ta frit 
was placed above the starting materials to act as a hlter 
during the centrifugation process. The Ta crucibles were 
subsequently sealed in quartz ampoules under ^1/3 atm 
Ar. 

The sealed ampoules were placed into a large box fur¬ 
nace and heated to 1190 °C over 12 h and held at that 
temperature for 24 h. The furnace was then cooled to 
800 °C over 390 h, after which the samples were removed 
and the flux was decanted using a centrifuge. Large, sin¬ 
gle crystal specimens (some on the order of ~1 g) of the 
Nd 2 Fei 4 B (2-14-1) phase (Fig. 1, inset) and smaller elon¬ 
gated prisms (~ 2 mm x 2 mm x 10 mm) of the NdFe 4 B 4 
(1-4-4) phase were extracted from the crucibles. The ul¬ 
timate sizes of the 2-14-1 crystals were not influenced by 
the compositions of the fluxes with the exception of the 
largest Ce concentration (in this case, large crystals of the 
Laves phase compound Ceo. 85 Ndo.i 5 Fe 2 formed rather 
than the 2-14-1 phase). The crystallographic faces of all 
phases present after decanting were well-defined. In Fig. 
1 we present XRD reflections resulting from an (0 0 1) 
face of the undoped Nd 2 Fei 4 B sample. 
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FIG. 1: XRD diffraction pattern from (0 0 1) face of undoped 
Nd 2 Fei 4 B crystal. The inset shows a size comparison image 
of crystals of the same composition. 
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Analyses of elemental compositions were performed 
using a Hitachi TM-3000 electron microscope equipped 
with a Bruker Quantax 70 energy dispersive (EDS) X-ray 
system (c/. Table I). For all compositions small crystals 
were ground into powder and used to produce XRD pat¬ 
terns using a PANalytical XPert Pro diffractometer (Cu 
Ka, 1.54056 A). The resultant peak reflections were ana¬ 
lyzed for lattice parameters using the LeBail fitting func¬ 
tion available with the onboard HighScore Plus software 
package. Single crystal diffraction was performed using 
a Rigaku single crystal X-ray diffractometer with Dectris 
Pilatus 200K detector (Mo Ka, 0.71073 A); the resul¬ 
tant diffraction patterns were analyzed using the onboard 
Rigaku software with absorption corrections performed 
using a spherical approximation. The structure was re¬ 
fined using SHELXL (Ref. [14]) with WinGX. Single 
crystal neutron diffraction was measured using the HB- 
3A four circle diffractometer at the High Flux Isotope 
Reactor (HEIR) at the Oak Ridge National Laboratory. 
The ~300 mg sample was measured at 4 K and 300 K 
by neutrons with a wavelength of 1.003 A from a bent 
Si-331 monochromator [15]. The structure was refined 
using the FULLPROF software package [16]. 

X-ray photoelectron spectroscopy (XPS) data were col¬ 
lected with a PHI 3056 XPS spectrometer with an A1 Ka 
source (1.4866 keV) in a cryo-pumped ultra-high vacuum 
chamber with a pressure of << I0“® Torr. Fresh surfaces 
of the crystal were created by grinding a single crystal in 
an Ar-filled glove box and transferring to the XPS cham¬ 
ber using a vacuum transfer system. Synchrotron data 
was taken at the Advanced Light Source (ALS) at the 
Lawrence Berkeley National Laboratory (Berkeley, USA) 
using Beamline 9.3.1. This beam-line is a bent magnet 
beam-line with an energy range of 2.3—5.2 keV. The min¬ 
imal spot size at the beam-line is 0.7 mm (v) x 1.0 mm 
(h) [17]. 

Magnetic properties were measured using 1) a Quan¬ 
tum Design Magnetic Property Measurement System 
(MPMS) and 2) the DC-extraction capability of the AC 
Measurement System (ACMS) option of a 14 T Quantum 
Design Physical Property Measurement System (PPMS). 
To measure the magnetization along the easy direction 
as a function of magnetic field 1) a c axis face was found 
on a single crystal using the X-ray diffractometer, 2) the 
crystal was polished into a parallelepiped with a long axis 
aligned with either the a or c direction (for hard and easy 
axis measurements, respectively) and 3) the crystal seg¬ 
ments were placed into the MPMS or PPMS under mag¬ 
netic fields of up to 13 T. Due to the tendency for the 
easy axis to align with the applied field direction, single 
crystal specimens were unable to quantitatively yield cor¬ 
rect hard axis magnetization measurements. To obtain 
these measurements, small crystals of each composition 
were ground into powder and placed into a gel cap with 
epoxy; the epoxy was allowed to set under an applied 
magnetic field so as to align the powders. Measurements 
were taken with H oriented perpendicular to the easy 
direction at 5 K, 300K, 350 K, and 400 K. To obtain 


TABLE II: Summary of 250 K single crystal X-ray diffraction 
data of sample (Ndo.78Ceo.22)2Fei4B; 

Space gronp P42/mnm; a = b = 8.8032(13)A; c = 
12.1880(20)A; a = /3 = 7 = 90° 


Atomic Coordinates 


Site 

X 

y 

z 

Uisa{AA 

B1 {4g) 

0.1239(8) 

0.1239(8) 

0 

0.0121(16) 

Pel (16fci) 

0.03732(7) 

0.35986(7) 

0.32411(5) 

0.0099(2) 

Fe2 (IGfcz) 

0.06698(7) 

0.27578(7) 

0.12753(5) 

0.099(2) 

Fe3 ( 8 ii) 

0.09802(7) 

0.09802(7) 

0.29599(7) 

0.0105(2) 

Fe4 ( 8 ^ 2 ) 

0.31758(7) 

0.31758(7) 

0.25407(7) 

0.0103(2) 

Fe5 (4e) 

0 

0 

0.11514(10) 

0.0098(3) 

Fe 6 (4c) 

0 

0.5 

0 

0.0105(3) 

Ndl/Cel ( 43 ) 0.22992(4) 

0.77008(4) 

0 

0.0108(2) 

Nd2/Ce2 (4/) 0.35743(4) 

0.35743(4) 

0 

0.0108(2) 

Reliability 





factors 

Ri 

WR 2 

Rint 

GOF 


0.0505 

0.1189 

0.0936 

0.882 


a correct value for the saturation magnetization a mea¬ 
surement of the same gel cap was taken with H oriented 
parallel to the easy direction; these data were then nor¬ 
malized to the values obtained from the single crystal 
easy axis measurements. The easy axis measurements of 
the aligned powder samples and the single crystal speci¬ 
mens were in good agreement. 

The M{T) properties of the 2-14-1 crystals were also 
measured using the Quantum Design MPMS. The mag¬ 
netization was measured as a function of temperature 
using an applied magnetic field of 10^ Oe; temperature 
ranges were 2 K < T < 750 K. Further confirmation 
of the Curie temperatures of the grown ferromagnetic 
crystals was obtained through thermal analysis using a 
Perkin Elmer Pyris Diamond Thermo-Gravimetric Ana¬ 
lyzer/ Differential Thermal Analyzer (TGA/DTA) with 
an applied magnetic field (< 0.1 T). 

The first principles calculations are carried out within 
density functional theory (DFT). We use the generalized 
gradient approximation [18] within the projector aug¬ 
mented wave method [19] as implemented in the Vienna 
ab initio simulation package [20, 21]. An energy cutoff 
of 400 eV and k space sampling on a 7 x 7 x 5 grid are 
employed. The Ce substitution on the Nd site is handled 
by the virtual crystal approximation (VGA). 

III. RESULTS AND DISCUSSION 
A. Crystal growth and Ce doping 

The as-grown crystals were roughly cubic in shape and 
displayed well-defined faceting (Fig. 1, inset). The sizes 
of the 2-14-1 crystals were typically ~ 2 x 2 x 2 mm^ with 
some reaching ~ 10 x 10 x 10 mm^ in volume. Except 
for the growth with a nominal flux composition of 18.1% 
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TABLE III: Summary of 250 K single crystal X-ray diffrac¬ 
tion data of sample (Ndo.78Ceo.22)2Fei4B; 

Space group P42/mnm; a = 6 = 8.8041(13)A; c = 
12.1633(17)A; a = ^ = 7 = 90° 


Atomic Coordinates 


Site 

X 

y 

z 

Uisa{°A^) 

B1 (4,) 

0.1221(14) 

0.1221(14) 

0 

0.008(3) 

Fel (16fci) 

0.03717(14) 

0.36003(14) 

0.32427(10) 

0.0093(4) 

Fe2 (I 6 /C 2 ) 

0.06714(14) 

0.27523(14) 

0.12785(11) 

0.0102(4) 

Fe3 ( 8 ji) 

0.09839(14) 

0.09839(14) 

0.29627(16) 

0.0108(4) 

Fe4 ( 8 ^ 2 ) 

0.31744(13) 

0.31744(13) 

0.25413(14) 

0.0102(4) 

Fe5 (4e) 

0 

0 

0.1155(2) 

0.0104(5) 

Fe 6 (4e) 

0 

0.5 

0 

0.0101(6) 

Ndl/Cel (4g) 

0.23015(7) 

0.76985(7) 

0 

0.0102(3) 

Nd2/Ce2 (4/) 

0.35750(7) 

0.35750(7) 

0 

0.0102(3) 

Reliability 





factors 

Ri 

WR 2 

Rint 

GOF 


0.0792 

0.1849 

0.0964 

1.196 


Nd and 32.9% Ce where a Laves phase compound formed 
rather than the 2-14-1 phase, all resulting single crystals 
were of similar size. The 1-4-4 phases grew as needles 
~ 1 X 1 X 10 mm^ and will be the subject of a forthcoming 
publication. 

Through EDS analysis, we were able to measure the 
relative concentrations of the rare earths (RE) and tran¬ 
sition metals (TM). The boron content could not be re¬ 
liably estimated from the EDS measurements. A variety 
of growths were used to maximize the Ce concentrations 
within the crystals. Results suggest that x = 0.38 is the 
maximum attainable by this synthesis route. 


B. Structural characterization 

The room temperature powder XRD patterns were 
taken of the 2-14-1 phases using the powder from ground 
crystals. LeBail fitting was used to extract the lattice pa¬ 
rameters, which are plotted in Fig. 2 as functions of Ce 
concentration. The values of both the a and c lattice pa¬ 
rameters decrease monotonically with increasing quanti¬ 
ties of Ce (c/. Fig. 2). At x = 0.38 in (Ndi_xCex) 2 Fei 4 B 
we see that the value of a is 0.14% smaller than that of 
the undoped sample. Similarly, for the same composi¬ 
tion, c is 0.22% smaller. No anomalies are present in the 
lattice parameters as a function of Ce content, in contrast 
with the sintered magnets of the same composition syn¬ 
thesized by Pathak et al. [4] where a two-phase region 
was observed for 0.15 < a: < 0.4 in (Ndi_xCex) 2 Fei 4 B. 
However, it should be noted that there is a wide spac¬ 
ing between our data points which could potentially belie 
small anomalies in the lattice parameters as a function of 
Ce content. Our single crystal specimens containing the 
two highest compositions of Ce {x = 0.22 and x = 0.38) 
fall within the two phase region suggested by the above 


authors, indicating that the range of solid solution in this 
system is sensitive to processing conditions such as tem¬ 
perature and solidification rate. 
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FIG. 2: Lattice parameters of the 2-14-1 phase plotted against 
Ce concentration. 

Single crystal X-ray diffraction on the sample with 
composition (Ndo. 78 Ceo. 22 ) 2 Fei 4 B was performed at 250 
K and 110 K; these temperatures bracket the spin re¬ 
orientation temperature of ~140 K in Nd 2 Fei 4 B [8, 22]. 
Structural refinements were made using SHELX [14] 
from 1068 and 966 unique reflections, respectively, and 
are presented in Tables 11 and 111. The space group 
matched that of previous reports and the powder diffrac¬ 
tion data, P4:2/mnm. The values from the single crystal 
data match fairly well with those obtained from powder 
diffraction. No discernible structural difference was seen 
between the data collected above the spin reorientation 
temperature and that collected below it. For the single 
crystal analyses no attempts were made to refine the RE 
position in terms of Nd/Ce occupancy as the scattering 
factors of the two RE elements were too similar. Later 
refinements using the occupancy results elucidated from 
the neutron diffraction data (below) yielded no change 
in the refinement. The structure derived from the 250 
K single crystal refinement is presented in Figs. 3a—3c 
and matches well with those published in the literature 
[6, 7, 23-25]. 

Single crystal neutron diffraction was performed on 
a ~300 mg sample of (Ndo. 78 Ceo. 22 ) 2 Fei 4 B at temper¬ 
atures of 300 K and 4 K. We collected ~630 reflections 
at each temperature. The magnetic signals were found at 
both 4 K and 300 K with a propagation vector of fc = 0, 
i.e., they coincide with the nuclear peaks. Therefore, in 
the data refinements, the nuclear and magnetic structure 
parameters were refined together. The results are pre¬ 
sented in Tables IV and V. 







5 



FIG. 3: 250 K structural refinement of sample 

(Ndo. 78 Ceo. 22 ) 2 Fei 4 B: a) isometric view showing the 

positions of the two different RE sites: dg and 4/ (the larger 
Ag site is denoted by a red label); b) view along a axis; c) 
view along c axis. 


Given the fact that the easy axis of this material is 
along the [0 0 1] direction, the magnetic vectors at 300 
K were constrained as such to reduce the total number 
of refined parameters. At 300 K the total refined mag¬ 
netic moment is 36.76±2.32 fig/i.u., in reasonable accord 
with the MPMS-determined value of 32.3 The Fe 

sites have moments ranging from 1.3(1)—3.2(2) ^b- This 
variability is likely caused by the different steric environ¬ 
ments experienced by each of the different Fe sites. With 
the singular exception of the 4e Fe site the size of the 
moments are correlated to the degree of Fe coordination 
around the Fe atoms. For example, two of the largest 
moments are seen on the 8 j 2 and 16fci sites which have 
Fe coordinations of 12 and 10, respectively. The smallest 
moment is observed at the Fe 4c site where the high¬ 
est RE coordination is present, with 4 of the 12 nearest 
neighbors being Nd/Ce. The two RE sites have moments 
of 2.372 and 1.605 fiB for 4(/ and 4/, respectively. The 
4^ site is, in terms of size, the larger of the two and 
is therefore predicted to be preferred for Ce, assuming 
the Ce is in the -1-3 oxidation state [9]. The Ce^’*' ion 
(115 pm diameter) is larger than the Ce'^+ ion (101 pm) 
and is therefore likely to prefer this larger site. When re¬ 
quiring that the thermal displacement parameters for the 
two different rare earth sites be equal, our 300K refine¬ 
ments presented in Table IV show that this preference 
is not absolute in that Ce site preference for the Ag site 
is 70%. If, however, we let the atomic displacement pa¬ 
rameters for each of the RE sites refine independently, 
the preference of Ce for the Ag site becomes even more 
absolute at 86 %. Though the degree of Ce preference for 
the Ag site and the atomic displacement parameters are 
highly correlated in these refinements, the fact that Ce 
does prefer the Ag site does not change. These observa¬ 
tions concur with the preliminary XPS data taken at the 
Advanced Light Source (ALS) at the Lawrence Berkeley 
National Laboratory whereby the XPS spectrum for the 


Ce 3d electrons most closely resembles that of Ce 203 , 
indicating that the majority of Ce is in the -1-3 oxidation 
state for this composition. Further experimentation is 
necessary to confirm this initial experimental result. 

For the 4 K refinement, we fixed the Ce occupancy 
to the values found at 300 K (70% on the RE Ag site 
and 30% on the RE 4/ site) and let the magnetic vec¬ 
tors vary in direction, as below ~140 K the easy axis is 
no longer the [0 0 1] direction. At 4 K, the total mag¬ 
netic moment for the sample is refined to be 49.1 ± 10.7 
/rsin reasonable agreement with the experimentally de¬ 
termined value of 34.9 /rs/f.u. for this composition (be¬ 
low) . The Fe sites have moments that refine in the range 
2.819—3.528 /is. Again, the moments associated with the 
Fe sites largely scale with their degrees of coordination 
by other Fe atoms. The refinement also shows that at 4 
K the moments are tilted by ~39.0°, consistent with the 
well-known phenomena of spin reorientation in the 2-14-1 
compounds at ^140 K [ 8 , 22] whereby the magnetic mo¬ 
ments orient themselves along the [3 3 5] direction rather 
than the [0 0 1 ] direction seen at room temperatures [26]. 
The angle between these two directions is 49.6°, vary¬ 
ing slightly from the experimentally determined value of 
39.0° from our refinement. However, to date no work 
has been done to quantify the effect of Ce substitution 
on this reorientation with respect to changes in the easy 
direction of magnetization. 


C. Spin re-orientation temperature 

Using single crystal neutron diffraction, the (0 0 6 ) re¬ 
flection of (Ndi_xCex) 2 Fei 4 B was measured upon warm¬ 
ing (Fig.4). Referring to Fig. 4, we see at 130.5 K a no¬ 
ticeable increase in the (0 0 6 ) peak intensity is evident, 
representing the spin reorientation transition whereby 
the magnetization easy axis shifts from the [0 0 1 ] direc¬ 
tion at high temperatures to something near the [3 3 5] 
direction at low temperatures. Previous investigations of 
undoped Nd 2 Fei 4 B by Abache and Oesterreicher [ 8 , 22] 
attributed this re-orientation to competing anisotropies 
of the Ag and 4/ RE sites. In their argument, the 4/ 
site (the smaller of the two) is more susceptible to pla¬ 
nar moment alignment. At temperatures above the spin 
re-orientation temperature (T^) the axial anisotropies of 
the RE and TM sites force the 4/ site to align axially 
in a ferromagnetic structure. However, below 140 K the 
preference of the 4/ site to align its moment in a planar 
orientation is much stronger compared to the competing 
influences of the Ee sites, resulting in spin re-orient at ion 
and a new easy axis for the system about 49° canted 
away from the c axis [ 8 , 22]. The spin re-orientation is 
therefore dependent on the 4/ site. The related com¬ 
pound Ce 2 Fei 4 B exhibits no spin re-orientation [ 8 ] im¬ 
plying that Ce has no inclination for planar orientation 
of its moment. Substituting Ce on the 4/ site is expected 
to affect the spin-reorientation. 

In Fig. 5a we plot the magnetization of undoped 
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FIG. 4: Changes in (006) reflection intensity as measured 
from 4-300 K using single crystal neutron diffraction on sam¬ 
ple with composition (Ndo. 78 Ceo. 22 ) 2 Fei/^B. 


Nd 2 Fei 4 B as a function of temperature and see increases 
in magnetization at 586 K and 143.5 K associated with 
Tc and Tg, respectively. In Fig. 5b we plot the evolu¬ 
tion of Tg as a function of Ce doping. From the neutron 
diffraction refinements we know that 18% of the RE 4/ 
site is substituted by Ce for the sample with the compo¬ 
sition (Ndo. 78 Ceo. 22 ) 2 Fei 4 B. Presumably similar levels of 
Ce are found at this site for the other Ce compositions 
investigated concomitant with the level of Ce introduced. 
We see from Fig. 5b that Tg is sensitive to the amount 
of Ce present, indicating that although Ce does prefer 
the ig site a not insignificant amount of Ce is still sub¬ 
stituting onto the 4/ site, assuming that the only way 
to modify the spin re-orientation temperature is to alter 
the occupancy of the 4/ site. This contrasts with the 
prediction of Alam et al. whereby Ce was predicted to 
sit on only the 4(; site [9]. 


D. Magnetic properties: Curie temperature, 
anisotropy field, and saturation magnetization 





o 
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X in (Ndj ,^Ce,j) 2 Fe, 4 B 


FIG. 5: a) Normalized Magnetization as a function of temper¬ 
ature for the undoped Nd 2 Fei 4 B sample taken using 1000 Oe 
showing both the Curie and spin re-orientation temperatures 
{Tc and Tg, respectively) and b) Evolution of spin reorienta¬ 
tion temperature as a function of Ce concentration. 


Relevant magnetic properties are summarized in Table 
IV. The Curie temperatures of the Ce-substituted 2-14- 
1 samples were measured via MPMS and DTA/TGA. A 
permanent magnet affixed to the exterior of the sample 
chamber were used to produce magnetic field gradient 
at the sample, so that a change in magnetic suscepti¬ 
bility produced a change in the apparent sample mass. 
The data from these analyses show a steady monotonic 
decrease in Tc with increasing Ce concentration. The 
undoped Nd 2 Fei 4 B yields a Tc of « 584 K. The compo¬ 
sition (Ndo. 62 Ceo. 38 ) 2 Fei 4 ^^B, on the other hand, shows a 


Tc of « 544 K, a total decrease of 1.1% Ce substituted. 
This result is less than the rate one would expect if a 
Vegard-like relation existed between the Nd 2 Fei 4 B and 
Ce 2 Fei 4 B (Tc « 430 K [8, 22]) parent compounds. 

In Figs. 6a and 6b we show the M{H) curves for the 
four Ce-substituted samples at 300 K and 400 K, respec¬ 
tively. With 38% Ce doping, the saturation magneti¬ 
zation, Ms, decreases 9.2% at 300 K (30.5 ^s/f.u. for 
(Ndo, 62 Ceo. 38 ) 2 Fei 4 B vs. 33.6 /rs/f.u. for the undoped 
Nd 2 Fei 4 B). Comparing these values to the 300 K sat- 
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uration magnetizations of the parent compounds (23.9 
/is/f.u. for Ce 2 Fei 4 B [27]) we find that the value of Mg 
seems to linearly scale with composition. The effects at 
higher temperatures are similar. At 400 K, for the same 
composition Mg is decreased by 7.4% from the undoped 
value. Concomitant with the decreases in Mg with in¬ 
creasing Ce concentration are decreases in Ha as deter¬ 
mined from extrapolating the hard axis magnetization 
curves in Figs. 6a and 6b to the saturation magnetiza¬ 
tion of the easy axis magnetization curves. At 300 K the 
value of Ha drops from an undoped value of 8.1 T to a 
value of 6.2 T for the highest achievable Ce concentra¬ 
tion. At 400 K Ha drops from 5.5 T to 4.7 T for crystals 
with large amounts of Ce present. 

IV. THEORETICAL CALCULATIONS 

The R 2 Fei 4 B class of materials (where R = Nd, Gd, Y 
etc.) has been the subject of much study [28-32]. With 
respect to Nd 2 Fei 4 B, first principles studies have been 
used to understand the effect of partial substitution of the 
Nd site by various rare earth elements (e.g., Y, La, Dy, Tb 
etc. [33, 34]) as well as the Fe site by Si, Ce, Sn etc. [35]. 
In this report we calculate the effects of Ce substitution 
for Nd on the magnetic properties of (Ndi_xCex) 2 Fei 4 B, 
focusing on the parameters relevant for a high perfor¬ 
mance magnetic material: the saturation magnetization. 
Curie temperature and magnetocrystalline anisotropy. 

A. Electronic and magnetic structnre of Nd2Fei4B 

Nd 2 Fei 4 B crystallizes in a tetragonal structure with 
four formula units per unit cell. In the atomic valence 
configuration the four electrons in the Nd-4/ orbitals and 
the six in the Fe-3ci orbitals make these ions strongly 
magnetic. The density of states (DOS) for Nd 2 Fei 4 B in 
the ferromagnetic state is shown in the blue curves in 
Fig. 7. Comparing the total DOS with partial DOS 
in the lower panels, we can see that the states near 
the Fermi level predominantly consist of Fe-d and Nd- 
/ states. There is an exchange splitting of ~2 eV in the 
Fe-d states with a nearly fully occupied spin-up chan¬ 
nel. The spin down channel has close to two electrons, 
suggesting a spin magnetic moment of ~3 sX the Fe 
sites, which is significantly above the value for bcc Fe 
(~2.2 /is/Fe). For comparison, the experimental neu¬ 
tron diffraction data collected at 4 K shows Fe having 
moments between 2.370 and 3.253 /rs, within the range 
of these extrema. 

At the Nd sites, the 4/ states are empty in the spin- 
up channel and are partially occupied in the spin-down 
channel, showing that in this calculation the Nd moments 
point in the opposite direction compared to the Fe mo¬ 
ments. The exchange splitting at the Nd-/ DOS is similar 
to that of Fe-d orbitals, even though the/ DOS is sub¬ 
stantially narrower. Counting the states in the Nd-/ / 




▼ (^^ 0 . 91 ^^. 089 ) 2 ^^ 14 ® 

■ (Nd„^gCeQ22)2Fsi4® 



b) —Nd,Fe,jB 
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(Nd,i.,sCe„ 22)2^^14® 

(Nd„_62Ce„.3g)2Fe|4B 



Magnetic Field, /UgH, T 


FIG. 6: M(H) curves for undoped and Ce-substituted 
Nd 2 Fei 4 B samples at: a) 300 K and b) 400 K. The insets 
display the anisotropy fields, Ha, as functions of doping. 


channel up to the Fermi level gives 3.1 electrons, which 
suggest a spin magnetic moment of -3.1 ^b at the Nd 
sites. From our calculations, we also find an orbital mo¬ 
ment of 2.8 at the Nd-/ states (opposite to the spin 
moment, as expected from Hund’s rules), while the or¬ 
bital moments at the Fe sites are ^0.06 /i^. The much 
larger orbital moment on the rare-earth atom suggests 



















Energy (eV) 

FIG. 7: The total and partial density of states for un¬ 
doped Nd-Fe-B (blue) and Ce doped Nd-Fe-B (red) from 
first-principles calculations in the ferromagnetic spin config¬ 
uration. The positive and negative DOS values correspond 
to spin-up and spin-down channels respectively. The DOS 
around Fermi level is dominated by Fe-d (middle panel) and 
Nd/Ce-/ (bottom panel) states. 


that the magneto-crystalline anisotropy is primarily a 
result of the rare-earth atom. Spin magnetic moments 
calculated by integrating charge density around the atom 
sites agree well with the above estimates. Integrating the 
spin-polarized charge density around the ions, we find av¬ 
erage moments of (in ^b) -3.4, 2.5 and -0.17 for Nd, Fe 
and B sites respectively. The Fe moment value is larger 
than that for pure Fe, contributing to the excellent mag¬ 
netic properties of this compound. For comparison, the 
experimental neutron diffraction data collected at 4 K 
shows Fe having net moments between 2.370 and 3.253 
^B, in good agreement with calculated value. Consis¬ 
tent with the experimental neutron diffraction data, we 
find that the Fe ions bonded with B atoms were found to 
have a reduced moment of 1.98 The net spin-|-orbital 
(J = S-l-L) magnetic moment of this system is 30.3 ^b 
per unit cell, in good agreement with both the neutron 
diffraction refinements and the magnetization data. 


B. Effects of Ce Substitution on Nd2Fei4B 

To understand the effect of Ce doping at the Nd sites 
on the properties of Nd 2 Fei 4 B, we replaced 10% of the 
Nd atoms by Ce using the virtual crystal approximation 
(VCA). The total and partial DOS for the doped system 
are presented as red curves in Fig. 7. As we can see. 


the only significant change in the DOS is at the RE-/ f 
states, which are shifted to lower energies by ~0.5 eV. 
Thus, we see that the electronic properties of the Nd-Fe- 
B magnets are not significantly affected by Ce doping, 
except a small reduction in the magnetic moments at the 
Nd/Ce site since Ce has two fewer electrons in its va¬ 
lence shells (4/^) than Nd (4/"^). Thus, the number of 
electrons in the 68 atom (8 rare earth atoms) unit cell is 
reduced by 1.6 electrons by 10% Ce substitution which 
lead to a subsequent reduction in the average magnetic 
moment at these sites to 3.2 fiB from 3.4 ^b- We find 
from these spin-polarized calculations that the Fe mo¬ 
ments are mostly unaffected by the substitution. 

An important consequence of the reduced exchange 
splitting is a reduction in Curie temperature, Tq- Gen¬ 
erally, for a ferromagnetic system characterized by large 
local moments such as the 2-14-1 materials, the Curie 
temperature is controlled by the difference in energy be¬ 
tween the ferromagnetic ground state and an essentially 
antiferromagnetic structure constructed so that most of 
the neighbors of Fe atoms are antiparallel. This in turn 
is a function of the exchange parameters connecting 
nearest neighbors i and j. In the simplest mean-field ap¬ 
proximation one may estimate the Curie temperature as 
a third of this energy difference, measured on a per Fe 
basis. 

We calculated the energy difference between the above- 
described ferromagnetic structure and an approximate 
antiferromagnetic structure obtained by flipping most of 
the nearest neighbor moments. The calculated energy 
difference is plotted in Fig. 8 as a function of the Ce 
concentrations. The x = 0 value of -2.71 eV per formula 
unit corresponds in the mean-field approximation to a 
Tc of 748 K, which is to be compared to the experimen¬ 
tal value of 585 K. The smaller experimental value sug¬ 
gests that effects beyond our simple mean-field approach 
may be relevant. A similar calculation overestimate was 
also found in Ref. [9]. Upon substituting Ce atoms via 
VCA, the magnitude of the energy difference decreases 
by about 4 percent for a 10 percent Ce substitution. This 
suggests that the mean field J decreases in Nd 2 Fei 4 B as 
X is increased which would in turn decrease the Curie 
temperature. This is consistent with our experimental 
results, which also find a decrease in Curie temperature 
with Ce alloying. 

Finally, we calculated the magnetic anisotropy energy 
(MAE) of (Ndi_xCex) 2 Eei 4 B for Ce concentrations x = 0 
and x = 0.1 to understand the effect of doping on this im¬ 
portant parameter. We used the experimental structure 
at room temperature (i.e., in the ferromagnetic state) to 
perform these calculations [36]. The total energy is then 
calculated in the presence of spin orbit coupling with the 
moments pointing along a and c directions. The MAE 
is defined as Ki = Ea — Ec- We find that for x = 0, 
Ki=12.3 meV/f.u. (^8 MJ/m^). This is somewhat 
larger than the experimental value of 4.3 MJ/m^, but 
we note that temperature effects are known to be signifi¬ 
cant in this compound [37]. We find that the anisotropy 
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Doping X 

FIG. 8: Variation of energy difference between ferromagnetic 
and antiferromagnetic structures in Nd 2 Fei 4 B as function of 
doping parameter x. The energy difference reduces upon Ce 
doping, suggesting that the exchange interaction J and con¬ 
sequently Tc will reduce as a function of Ce concentration. 

energy remains more or less same upon doping; we get 
Ki=13.2 meV/f.u. when x = 0.1. We do not consider 
the increase to be significant, but expect that larger lev¬ 
els of Ce alloying would tend to show a reduction in the 
value of the MAE. In short, we find that that the high 
magnetic anisotropy energy of the Nd-Fe-B system is not 
significantly affected by the Ce doping, which confirms 
experimental observations. 

V. SUMMARY 

Large single crystals of Nd 2 Fei 4 B were grown from 
Nd-Fe flux. These crystals were a maximum of ^0.5 


g in mass and displayed sharp faceting. Cerium sub¬ 
stitution was attempted for Nd; the maximum degree 
of substitution we were able to achieve was the com¬ 
position (Ndo, 62 Ceo, 38 ) 2 Fei 4 B. We found that Ce and 
Nd formed a continuous solid solution in contrast to the 
melt-spun ribbons of similar compositions fabricated by 
Pathak et al. [4] where a two phase region was observed 
for 0.15 < X < OA in (Ndi_xCex) 2 Fei 4 B. Single crystal 
neutron diffraction showed that the Ce has a slight pref¬ 
erence for the larger RE Ag site, favoring this location by 
a factor of ~7:3 over the RE 4/ site. 

The presence of Ce serves to deleteriously affect mag¬ 
netic properties such as Tc, Mg, and Ha. The values 
of these characteristics largely fall in line with empirical 
predictions based on the application of simple alloying 
rules between the two parent compounds (Nd 2 Fei 4 B and 
Ce2Fei4B). 

We have also used X-ray photoelectron spectroscopy 
and synchrotron emissions to determine the oxida¬ 
tion states and binding energies of all species in Ce- 
substituted Nd 2 Fei 4 B; preliminary data shows that the 
XPS spectrum for the Ce M electrons most closely re¬ 
sembles that of Ce 203 . 

With the help of first-principles calculations, we stud¬ 
ied the electronic and magnetic properties of Ce doped 
Nd-Fe-B magnets. We find that for up to a; = 0.1 dop¬ 
ing of Ce, the electronic structure around the Fermi level 
is not significantly affected. The exchange splitting and 
the magnetic moments at the Nd sites are diminished 
as a consequence of Ce doping, which in turn reduces 
the Curie temperature. Finally, we find that the high 
magnetic anisotropy energy of the Nd-Fe-B system is not 
significantly altered by the Ce doping, confirming the ex¬ 
perimental observations. 


[1] Critical Materials Strategy. Technical report, United 
States Department of Energy, December 2011. 

[2] M. Sagawa, S. Fujimura, H. Yamamoto, Y. Matsuura, 
and K. Hiraga. Permanent Magnet Materials Based 
on the Rare Earth-Iron-Boron Tetragonal Compounds. 
IEEE Transactions on Magnetics, MAG-20:1584-9, 1984. 

[3] J. J. Croat, J. F. Herbst, R. W. Lee, and F. E. Pinker¬ 
ton. Praseodymium-Iron- and Neodymium-Iron-Based 
Materials: A New Class of High-Performance Permanent 
Magnets. Journal of Applied Physics, 55:2078-82, 1984. 

[4] A. K. Pathak, M. Khan, . A. Gschneidner, Jr., R. W. Mc- 
Callum, L. Zhou, K. Sun, K. W. Dennis, C. Zhou, F. E. 
Pinkerton, M. J. Kramer, and V. K. Pecharsky. Cerium: 
An Unlikely Replacement of Dysprosium in High Perfor¬ 
mance Nd-Fe-B Permanent Magnets. Advanced Materi¬ 
als, page Ahead of Print, 2015. 

[5] D. Givord, H. S. Li, and J. M. Moreau. Magnetic 
Properties and Crystal Structure of Neodymium-Iron- 


Boron (Nd 2 Fei 4 B). Solid State Communications, 50:497- 
9, 1984. 

[6] J. F. Herbst, J. J. Croat, F. E. Pinkerton, and 
W. B. Yelon. Relationships Between Crystal Struc¬ 
ture and Magnetic Properties in Neodymium-Iron-Boron 
(Nd 2 Fei 4 B). Physical Review B, 29:4176-8, 1984. 

[7] J. F. Herbst. Rare Earth(R) Iron Boron (R 2 Fei 4 B) Ma¬ 
terials: Intrinsic Properties and Technological Aspects. 
Reviews of Modern Physics, 63:819-98, 1991. 

[8] C. Abache and J. Oesterreicher. Magnetic Anisotropies 
and Spin Reorientations of R 2 Fei 4 B-Type Compounds. 
Journal of Applied Physics, 60:3671-9, 1986. 

[9] A. Alam, M. Khan, R. W. McCallum, and D. D. John¬ 
son. Site-Preference and Valency for Rare-Earth Sites 
in (R-Ce) 2 Fei 4 B Magnets. Applied Physics Letters, 
102:042402/1-042402/4, 2013. 

[10] P. C. Canfield and 1. R. Fisher. High-Temperature So¬ 
lution Growth of Intermetallic Single Crystals and Qua- 




10 


TABLE IV: Summary of 300 K single crystal neutron diffraction data of sample (Ndo. 78 Ceo. 22 ) 2 Fei 4 B 


Atomic Coordinates 


Lattice Constants 

Site 

X 

y 

z 

Uisa{A^) 

Occ. 
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B1 {4g) 

0.12430(33) 

0.12430(33) 

0 

0 . 0110 ( 8 ) 

1 

c = 12.18801 

Fel (16fci) 
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0.35957(13) 
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0.0081(3) 

1 

a = /3 = 7 = 90° 
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0.0081(3) 

1 
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1 

PJ^franm 
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Ce2 (4/) 
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0.13 



Magnetic vectors 





r ihs) 

<!> 

e 



Total Moment 

Fel (16fci) 
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0 

0 



36.76 ±2.32 pg 

Fe2 ( 16 fc 2 ) 

2.247(135) 

0 

0 




Fe3 ( 8 ji) 

1.614(169) 

0 

0 




Fe4 ( 832 ) 
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0 

0 




Fe5 (4e) 

3.247(159) 

0 

0 




Fe6 (4c) 

1.305(137) 

0 

0 




Ndl ( 43 ) 

2.372(119) 

0 

0 




Cel ( 43 ) 

2.372(119) 

0 

0 




Nd2 (4/) 

1.605(118) 

0 

0 




Ce2 (4/) 

1.605(118) 

0 

0 



Reliability factors 

Rf2 

Wf2 

Rf 
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0.0461 

4.06 
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TABLE V: Summary of 4 K single crystal neutron diffraction data of sample (Ndo.78Ceo.22)2Fei4B 


Atomic Coordinates 


Lattice Constants 

Site 

X 

V 

z 

C/i«o(A2) 

Occ. 

a = b = 8 . 8 OOOA 

B1 (4p) 

0.12377(41) 

0.12377(41) 

0 

0.0051(11) 

1 

c = I 2 .I 88 OA 

Fel (16fci) 

0.03736(16) 

0.36009(16) 

0.32426(11) 

0.0018(4) 

1 

a = ^3 = 7 = 90° 

Fe2 (16^2) 

0.06698(17) 

0.27540(17) 

0.12759(11) 

0.0018A) 

1 


Fe3 (8/1) 

0.09795(17) 

0.09795(17) 

0.29655(17) 

0.0018A) 

1 

PA2/mnm 

Fe4 (8/2) 

0.31775(17) 

0.31775(17) 

0.25415(16) 

0.0018A) 

1 


Fe5 (4e) 

0 

0 

0.11597(24) 

0.0018A) 

1 


Fe6 (4c) 

0 


0 

0.0018A) 

1 


Ndl (4p) 

0.22988(28) 

0.77012(28) 

0 

0.0009(5) 

0.69 


Cel {4g) 

0.22988(28) 

0.77012(28) 

0 

0.0009(5) 

0.31 


Nd2 (4/) 

0.35713(27) 

0.35713(27) 

0 

0.0009(5) 

0.87 


Ce2 (4/) 

0.35713(27) 

0.35713(27) 

0 

0.0009(5) 

0.13 



Magnetic vectors 





r (Ms) 

<t> 

e 



Total Moment 

Fel (16fci) 

3.288(76) 

226.09(569) 

39.00(180) 



36.76 ±2.32 pB 

Fe2 (16^2) 

2.878(76) 

199.04(510) 

39.00(180) 




Fe3 (8/1) 

3.378(115) 

156.65(527) 

39.00(180) 




Fe4 (8/2) 

3.528(92) 

191.39(471) 

39.00(180) 




Fe5 (4e) 

3.302(103) 

240.22(706) 

39.00(180) 




Fe6 (4c) 

2.819(131) 

210.84(711) 

39.00(180) 




Ndl (Ag) 

2.921(70) 

255.63(732) 

39.00(180) 




Cel (Ag) 

2.921(70) 

255.63(732) 

39.00(180) 




Nd2 (4/) 

3.644(96) 

185.29(536) 

39.00(180) 




Ce2 (4/) 

3.644(96) 

185.29(536) 

39.00(180) 



Reliability factors 

Rf2 

Wf2 

Rf 





0.0459 

0.0948 

0.0648 

0.0738 
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